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Abstract 

The compound LiAl^T^-yC^ undergoes a metal-to-insulator transition for y c ~ 0.33. It is 
known that disorder alone is insufficient to explain this transition; e.g., a quantum site percolation 
model predicts y c ~ 0.8. We have included (Hubbard) electronic interactions into a model of this 
compound, using a real-space Hartree-Fock approach that acheives self consistency at every site, 
and have found that for a Hubbard energy equal to 1.5 times the non-interacting bandwidth one 
obtains y c ~ 0.3. Further, with increasing Hubbard energy we find an Altshuler-Aronov suppression 
of the density of states, 5N(e) ~ yje — e_p|, that reduces the density of states at the Fermi energy 
to zero at the critical Hubbard interaction. Using this ratio of correlation to hopping energy one 
is led to a prediction of the near-neighbour superexchange (J/t ~ 1/3) which is similar to that for 
the cuprate superconductors. 
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The range of interesting phenomena of transition metal systems includes the anomalous 
normal state and high-temperature superconductivity of the quasi-2d cuprates Q], as well 
as the colossal magnetoresistance of the manganites and related systems j^j. The ubiquitous 
physics believed to be responsible for these novel behaviours is strong electron-electron 
interactions, a consequence of which is the inapplicability of Landau's theory of fermi liquids. 
One interesting class of related materials are the spinels. There are over 300 transition metal 
spinels, out of which only four are superconducting. Further, of these four only one is an 
oxide, and that oxide, LiTi 2 4 , has the highest superconducting transition temperature, 
T c 0: CuRh 2 Se 4 (T c = 3.49 K), CuV 2 S 4 (T c = 4.45 K), CuRh 2 S 4 (T c = 4.8 K), and 
LiTi 2 4 (T c = 11.3 K). 

LiTi 2 4 is a 1 / 4-filled d ' 5 system in which the electronic conduction occurs via direct 
d — d hopping on the Ti sites, owing to the orientation of the low-lying t 2 g orbitalsj^O]. The 
Ti sublattice corresponds to a corner-sharing tetrahedral lattice (CSTL), which is a fu 
frustrated three dimensional (3d) structure. It has been suggested by Bednorz and Miiller 
that this system is moderately correlated electronic system, and that the superconductivity 
may be driven by the electronic interactions amongst the d electrons. However, other reasons 
behind our interest in this system include: (i) A metal-to-insulator transition generated by 
both disorder and a reduced electronic density caused by chemically substituting Li, Al or 
Cr for Ti this transition is the focus of this paper, (ii) If strong electronic correlations 
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are present then the magnetic properties of this material could be interesting, corresponding 

in) The isostructural d 1 ' 5 LiV 2 4 is believed to be 



to a 1 / 4-filled fully frustrated lattice. 

n 

a d-electron heavy fermion compound |9j, and understanding the simpler (1/2 an itinerant 
electron per site in LiTi 2 4 vs. the 1.5 electron per site (1 local moment plus 1/2 an itinerant 
electron per site)) LiV 2 4 would be highly beneficial. 

In the tight-binding approximation, a reasonable Hamiltonian from which to begin a 
study of disordered LiAl ?; Ti 2 _ ?/ 4 is given by 

H = ^ EjU^a - t ^2(4 a c ja + h.c.) + U '^n^Uii (1) 

i,cr i 

where % denotes a Ti site on a CSTL, and Cj CT the annihilation operator for an electron at 
site i with spin a, and rij jCr = c ia Ci a . (For reference below, note that the noninteracting 
bandwidth is 8t.) 

As mentioned above, there are a variety of different chemical dopings that lead to a 



metal-to-insulator transition (MIT), but here we focus on LiAl^T^-yC^, which undergoes 
a MIT for y ~ 0.33 gj. In this system the Al 3+ ions substitute onto the Ti sublattice, and 
thus in a first approximation the on-site energies are chosen at random according to the 
distribution function 

P{Si) = (1 " |) 6(£i - £ T i) + | S{6i ~ £ A l) • (2) 

Then, in the expected limit of Eai — ^Ti ^ 8t, if one ignores the presence of electron-electron 
interactions one sees that this system is an excellent representation of quantum site percola- 
tion (QSP). Previously, three of us have analyzed such a QSP model for noninteracting 
electrons (U = 0) on a CSTL, and (numerically) exactly solved the disordered electron 
problem (energies and eigenfunctions) for various large lattice sizes, thus determining the 
Fermi energy, Ef(y), and mobility edge, E c (y), as a function of Al doping concentrations. 
We found that at y = y c ~ 0.8 the Fermi energy and mobility edges crosses, and in such 
a disordered, noninteracting model this would correspond to the predicted MIT. The large 
disagreement between the experimental value of y c ~ 0.33 and our prediction (y c ~ 0.8) 
highlights that important physics has been omitted in such an analysis [ill ]. 

Indirectly, the above result supports the conjecture that strong electronic correlations are 
present and play an important role in the physics of LiTi 2 04. Of course, to provide substance 
to this idea we need to examine the full disordered and interacting electron problem, and 
to this end we have analyzed the behaviour of such a disordered-Hubbard Hamiltonian, viz. 
we have studied Eq. {1} with the Ti conduction path limited by the full QSP model men- 
tioned above, now including the Hubbard interaction term. This is a relatively complicated 
model which contains the interplay between electronic correlations and disorder produced 
by quantum-site-percolation, and in the following we describe our results of a comprehen- 
sive examination of this problem in a real-space self-consistent Hartree-Fock approximation. 



Recent results |l2j for a disordered Hubbard model on a two-dimensional square lattice have 
lead to the interesting prediction of a novel metallic phase; here we use the same formalism, 
but now for a 3d CSTL. 

In such a real-space self-consistent Hartree-Fock formulation one replaces the Hubbard 
interaction term as follows: 

U^2 n iT n a ?iU^2 (Ki)^ + (n^n^ - (n^n^)) , (3) 



where the primed summation indicates that only Ti sites (from both the maximally con- 
nected, and possibly isolated clusters) are included. Then the interacting Hamiltonian is 
simplified to 

K//= ~t ^'clc^-Uj^in^in^ (4) 
+ Uj^ [(n^)c\ lCii + (n^cl^} 

i 

where 

M= \<?\M\ 2 , (5) 

a 

and 

T-(-eff\^aa) = £aa\^Paa) (6) 

Here, |i) is a single orbital Wannier function at site % from the Ti sublattice, \ip aa ) are the 
energy eigenstates of the effective Hamiltonian, (ni a ) is the average density of electrons with 
spin a at site i, and Ep is the Fermi energy |l3l |. 

In order to find self-consistent solutions for the above Hamiltonian, we start with the 
noninteracting but disordered tight-binding model on a CSTL, and obtain the expectation 
values of the density of electrons on each site using Eq. (J3J) - this is the distribution of the 
density of electrons for a non-interacting QSP disordered system for a specific filling factor 
of the LiAl ?/ Ti2-y04 system. We then iterate to convergence at every site as U is increased 
from zero. Note that in the iteration sequence we have added a small random fluctuation 
in the on-site densities for each site and each spin, to allow for the system to proceed to a 
non-paramagnetic state if it so chooses. 

This procedure has been applied to CSTL, with random QSP determined by y, for system 
sizes of 5488, 8192, 11664 and 16000 lattice sites with periodic boundary conditions. For 
each system size the spectrum of eigenvalues and the distribution of the electron density 
throughout the lattice were calculated for each realization of disorder, for doping concen- 
trations of y = 0.2,0.3,0.4,0.6,0.7, and 0.75, and for Hubbard interaction strengths of 
U/t — 2, 4, 6, 8, 10, and 12. Then, these calculations were repeated for different complexions 
of disorder for each set of values of y and U/t. 

Density of states results are shown in Fig.^for a doping concentration of y = 0.3 and for 
three different Hubbard interactions of U/t = 4,8, 12, for systems with 16000 lattice sites. 
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FIG. 1: Density of states for a single realization of disorder in the Hubbard+QSP model for the 
LiAl y Ti2- J/ 04 system with y = 0.3 and U/t = 4,8, 12 for a system with 16000 lattice sites. The 
suppression that is associated with the Altshuler-Aronov square-root singularity is evident, as is 
the complete suppression (zero density of states at the Fermi level) at the critical (U c /t,y c ). 



The interesting feature in these plots is the suppression of the density of states at the Fermi 
level with increasing U. That is, as a function of increasing Hubbard energy, we find that 
a progressively larger change of the density of states consistent with the predictions of the 
theory of Altshuler and Aronov 3], namely a square- root suppression near the Fermi level 
of the form SN(e) ~ yje — ep\. In fact, we find a complete suppression that first appears at 
the critical Hubbard interaction (see below) associated with the metal-to-insulator transition. 

n 

These and our previous density of states data|H| are consistent with the hypothesis that 
only in an interacting system should one find such effects. Lastly, we note that recent 
experimental investigations of other transition metal oxides have seen precisely this 
form of the density of states, on the metallic side, as the metal-to-insulator transition is 
approached. 

In order to characterize the metallic or insulating behaviour of such systems, we exam- 
ined the eigenstates of a small energy bin (AE/t = 0.1) located symmetrically about the 
Fermi energy, and calculated the average inverse participation ratio (IPR), averaged over a 
sufficient number of realizations to obtain converged data. A collection of our results for the 
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FIG. 2: Plots of the IPRs vs. the inverse of the system size for the Hubbard+QSP model for 
the LiAl y Ti2_j / 04 system with y = 0.2, 0.3 and 0.4 for different strengths of the Hubbard on-site 
interaction. Note that the vertical scale in all of these plots is identical. 



IPRs vs. the inverse of the system size is shown in Fig. |21 From the linear fits of the IPRs 
vs. the inverse of the system size, the values of the IPRs in the thermodynamic limit can be 
extracted. We also determined the 98% confidence level (2x the standard deviation) error 
bars of the intercepts. As is well known, if the extrapolated value of the IPR has a finite 
intercept then the eigenstates corresponding to that specific system, in that energy range, 
are localized, and in the subsequent discussion we use the criterion that the intercept must 
be at least twice the standard deviation above zero before we classify (with a 98% confidence 
level) that those energy eigenstates are of a localized nature. For example, for y = 0.3 for 
U/t = 2 to 10 the states near the Fermi level are extended, whereas for U/t = 12 these 
same states are localized. Consequently, as a function of increasing Hubbard interaction we 
identify the MIT as occurring at U c /t ~ 11 — 12 for y = 0.3. 

Repeating this sequence of calculations for the above-mentioned y and U/t we identified 
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a phase diagram of the (T = 0) metal-to-insulator transition of this system with respect 
to disorder and interaction, and our results are depicted in Fig. 3. We immediately see 
that if strong electronic correlations are indeed an important aspect of the physics associ- 
ated with this transition, that is for a y e c x P l ~ 0.33, we require a U/t ~ 11 — 12. Noting 
that for the noninteracting electrons on this lattice one has a bandwidth of 8t, this implies 
that LiTi204 is a moderately correlated three-dimensional electronic system. This is con- 
sistent with conclusions drawn, albeit indirectly, from recent low-temperature specific heat 
measurements in this material. 

In fact, this value is not far from experimental estimates of U/t for the t2 9 d orbitals 
of related Ti-based transition metal oxides. That is, the estimate for the Hubbard on-site 
repulsion for the t2 9 orbitals of Ti atoms in the perovskite LaTiC>3 is U t2g ~ 3.1eV jl7|. 
Resonant soft-x-ray emission spectroscopy on perovskite-type Ti compound La x Sri_ x Ti03 
provides an estimate of Ud-d ~ 4.0 — AAeV [la . Q], with the most recent estimate being 



U/t ~ 4.(W(2003)|12|. Now compare these estimates to our predicted value of U for 
the LiALTi9_ u 04 system: taking into account the estimate for the transfer integral of 

nn 

t ~ 0.33eV for the t2 9 band in this system based on the LDA calculations |4J, la], using 
U c /t = 11 — 12 (for y = 0.3) our estimate of U corresponds to 3.7-4.0 eV. Clearly, our 
estimate is not dissimilar to the experimentally observed values for these related Ti-based 
oxides. 

We note that this energy leads to a provocative comparison between the near-neighbour 
superexchange (J) between moments in this system vs. those in the high-temperature 
superconducting cuprates. That is, using J = 4t 2 /U for this one-band system we find 
that in LiTi 2 04 J/t m 1/3, which is similar to the estimates for the cuprates. So, this 
similarity also lends support to the hypothesis of Bednorz and Miiller |(| of the potential 
relation of the pairing in this system to the cuprates. 

As further corroboration of our phase diagram results, and to gain a better understanding 
of how the transition is connected to the electronic and magnetic properties of this system, 
we have examined the (spin-resolved) charge and magnetic densities as a function of disorder 
and Hubbard energy. Our results for the variation of the number of electrons per site are 
effectively independent of spin, and are shown in Fig. 4. We see that for y = 0.3, for U/t = 
to 8, essentially no change takes place. However, as U/t is increased from 8 to 10 to 12 this 
quantity undergoes substantial change. That is, when the Hubbard energy approaches the 
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FIG. 3: Phase diagram for the metal-insulator transition in the Hubbard+QSP model represen- 
tative of the LiAl J/ Ti2_j / 04 system. Filled circles are estimates from the thermodynamic limit 
extrapolation of the IPRs, discussed in the text. 

critical value (for this y) due to the proliferation of localized states arising from the suppres- 
sion of the density of states at the Fermi level, one finds a far more inhomogeneous system. 
Also, our numerical results for the local magnetizations correspond to the appearance of 
(short-ranged) antiferromagnetic correlations at the metal-to-insulator transition. That is, 
as U/t is increased towards U c /t(y) this distribution gradually changes from a rather narrow 
peak at zero (paramagnetism), to a distribution dominated by two broad peaks at ±mo, 
with m ~ 0.4, corresponding to a locally antiferromagnetic arrangement but with many 
sites having largely paramagnetic character. 

Summarizing, we have used a real-space self-consistent Hartree-Fock formulation to ex- 
amine the metal-to-insulator transition in disordered LiAl J/ Ti 2 _j / 04; this work (i) treats the 
disorder exactly (that is, diagonalizing the full Hamiltonian matrix for any particular com- 
plexion of disorder), and (ii) involves solving for the self-consistent solutions at every site. 
We have found that to obtain agreement with the experimentally observed concentration of 
Al impurities at which the metal-to-insulator transition occurs, a Hubbard energy somewhat 
larger than the noninteracting bandwidth is required, consistent with recent experiments on 
this system 13], and that the resulting density of states as the transition is approached is 
similar to that found in experiments on related materials jl^ ]. 

We thank Nandini Trivedi and Barry Wells for helpful comments. This work was partially 



8 



0.2 0.4 0.6 0.8 

Average # of electrons per site 



FIG. 4: The distribution of single-site charge densities for y = 0.3 as a function of increasing 
Hubbard energy U/t. One sees a profound change in the homogeneity of the charge density as U jt 
approaches the critical value at which the metal-to-insulator transition occurs. 
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